Huge Positive Magnetoresistance in Antiferromagnetic Double Perovskite Metals 
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Metals with large positive magnetoresistance are rare. We demonstrate that antiferromagnetic 
metallic states, as have been predicted for the double perovskites, are excellent candidates for huge 
positive magnetoresistance. An applied field suppresses long range antiferromagnetic order leading 
to a state with short range antiferromagnetic correlations that generate strong electronic scattering. 
The field induced resistance ratio can be more than tenfold, at moderate field, in a structurally 
ordered system, and continues to be almost twofold even in systems with ~ 25% antisite disorder. 
Although our explicit demonstration is in the context of a two dimensional spin-fermion model of the 
double perovskites, the mechanism we uncover is far more general, complementary to the colossal 
negative magnetoresistance process, and would operate in other local moment metals that show a 
field driven suppression of non-ferromagnetic order. 



There has been intense focus over the last two decades 
on magnetic materials which display large negative mag- 
netoresistance (MR) [THS]- In these systems, typically, an 
applied magnetic field reduces the spin disorder leading 
to a suppression of the resistivity. The field may even 
drive an insulator-metal transition leading to 'colossal' 
MR ;lj. Large positive MR is rarer, and seems counter- 
intuitive since an applied field should reduce magnetic 
disorder and enhance conductivity. We illustrate a sit- 
uation in double perovskite [4 metals where an applied 
field can lead to enormous positive MR. The underlying 
principle suggests that local moment antiferromagnetic 
metals [MH], at strong coupling, should in general be 
good candidates for such unusual field response. 

The double perovskites (DP), A2BB'06, have been ex- 
plored 0] mainly as ferromagnetic metals (FMM) or anti- 
ferromagnetic insulators (AFMI). The prominent exam- 
ple among the former is Sr2FeMo06 (SFMO) fS^, while 
a typical example of the later is Sr2FeW06. The FMM 
have moderate negative MR [5Mllj. Like in other corre- 
lated oxides [12], the magnetic order in the ground state 
is expected to be sensitive to electron doping and it has 
been suggested P^HT5] that the FMM can give way to 
an AFM metal (AFMM) on increasing electron density. 
The AFM order is driven by electron delocalisation and, 
typically, has lower spatial symmetry than the parent 
structure. The conduction path in the AFM background 
is 'low dimensional' and easily disrupted. 

There is ongoing effort [111 [TTj , but no clear indication 
yet of the occurence of an AFMM on electron doping 
the FMM double perovskite. The problems are twofold: 
(a) Antisite defects: in a material like SFMO, substitut- 
ing La for Sr to achieve a higher electron density in the 
Fe-Mo subsytem tends to make the Fe and Mo ionic sizes 
more similar (due to resultant valence change), increases 
the likelihood of antisite disorder (ASD), and suppresses 
magnetic order, (b) Detection: even if an AFM state is 
achieved, confirming the magnetic order is not possible 
without neutron scattering. The zero field resistivity is 
unfortunately quite similar T^ to that of the FMM. 



We discovered a remarkably simple indicator of an 
AFMM state. Using a real space Monte Carlo method 
on a two dimensional (2D) model of DP's we observe 
that: (i) For temperatures below the zero field Tc of the 
AFMM, a modest magnetic field can enhance the resis- 
tivity more than tenfold, (ii) The magnetoresistance is 
suppressed by structural defects (antisite disorder) but 
is still almost twofold for ~ 25% mislocation of B sites, 
(iii) The enhancement in resistivity is related to the pres- 
ence of short range correlated AFM domains that survive 
the field suppression of global AFM order and lead to 
enhanced electron scattering. The mechanism suggests a 
generalisation to 3D, where direct simulations arc diffi- 
cult, and to other local moment AFMM systems where 
an applied field can suppress long range AFM order. 

We first describe the structurally ordered DP model 
and will generalise it to the disordered case later. 
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f is the electron operator on the magnetic, B, site and m 
is the operator on the non-magnetic, B', site, ts and 
tB' are onsite energies, at the B and B' sites respec- 
tively. The B and B' alternate along each cartesian axes. 
A = e_B — eB' is a 'charge transfer' energy, t is the hop- 
ping between nearest neighbour (NN) B and B' ions on 
the square lattice and we will set t = 1 as the reference 
scale. Si is the core B spin on the site Ri. We assume 
S = I Si I — 1 and absorb the magnitude of the spin in 
the Hund's coupling J at the B site. We use J/t ^ 1. 
When the 'up spin' core levels are fully filled, as for Fe in 
SFMO, the conduction electron is forced to be antipar- 
allel to the core spin. We have used J > to model 
this situation. We have set the effective level difference 
Ae// = (es — J/2) — tB' = 0. The chemical potential 
/i is chosen so that the electron density is in the window 
for A type (stripe like) order of core spins. N is the total 
electron count and h is the applied magnetic field in the z 



direction. We have ignored next neighbour hopping and 
orbital degeneracy in the present model. We also focus 
on the 2D case here, where simulation and visualisation 
is easier, and comment on the 3D case at the end. 

The model has a FMM ground state at low electron 
density (TSHU]- This is followed by a phase with stripe- 
like order - with FM B lines coupled AFM in the trans- 
verse direction. We call this the 'A type' phase. This is 
followed by a more traditional AFM phase where an up 
spin B ion, say, is surrounded by four down spin B ions, 
and vice versa (the 'G type' phase). We focus here on 
the A type phase since it has a simple 3D counterpart 
and occurs at physically accessible electron density 

We solve the spin-fermion problem via real space 
Monte Carlo using the traveling cluster method [H] . This 
allows access to system size 32 x 32. We have used field 
cooling (FC) as well as zero field cooling (ZFC) protocols. 
For ZFC the system is cooled to the target temperature 
at h/t = and then a field is applied. We calculate the 
resistivity and the magnetic structure factor peaks and 
also keep track of spatial configurations of spins. The 
optical conductivity is calculated via the the Kubo for- 
mulation, computing the matrix elements of the current 
operator. The 'dc conductivity' is the low frequency av- 
erage, adc = (1/Aw) /q a{uj)dLL>. This is averaged over 
thermal configurations and disorder, as appropriate. Our 
'dc resistivity' p is the inverse of this adc- Aw = 0.05i. 

The A type pattern has two possible orientations of 
the stripes, either from bottom left to top right, or from 
bottom right to top left. These are the two 'diagonals' 
in 2D. The first corresponds to peaks in the structure 
factor D{q) at {Qai,Qa2}, and the second to peaks at 
{QA3,QAi}- For reference, Q^i = {7r/2,7r/2}, Qa2 = 
{37r/2,37r/2}, Qa3 = {7t/2,3tt/2}, Qa4 = {3^2,^/2}. 
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FIG. 1: Colour online: The resistivity p{T) in the absence of 
antisite disorder for cooling in different applied fields. For 
temperatures below the zero field transition, T^ , p{T) in- 
creases on applying a field, and for T > T^ p{T) decreases 
on applying a field. The ratio p{T,h)/p(T,0) can be very 
large as T — > in this field cooling situation. 
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FIG. 2: Colour online: The field dependence of the resistivity 
at different temperatures following a zero field cooling (ZFC) 
protocol. For T = T^ (just above Tc) the zero field resistivity 
is already large and p decreases slightly with h due to sup- 
pression of spin disorder. For T — T~ (just below Tc) and at 
Tc/2 there is a sharp increase in resistivity at h/t ~ 0.01, with 
a peak around h/t ~ 0.02 — 0.03 and a fall thereafter. This 
is consistent with the trends seen in Fig.l. At T = Tc/4 this 
ZFC scheme does not manage to create competing magnetic 
structures for h/t ~ 0.05, possibly due to metastability of the 
parent AFM pattern. The field induced switching is therefore 
easiest achieved between Tc/2 and Tc- 



The ferromagnetic peaks in S'(q) are at Q^^i = {0, 0} and 
Qf2 = {7r,7r}. The ordered configurations lead to peaks 
at two wavevectors since the model has both magnetic 
and non-magnetic sites and our wavevectors are defined 
on the overall B-B' lattice. 

Let us first examine the FC results in resistivity, Fig.l. 
Cooling at h/t = leads to a sharp drop in resistivity 
at T = r° ~ 0.032, where T° is the zero field transition 
temperature, and p{T) — > as T — > 0. Cooling at h/t = 
0.01 leads to a small suppression in Tj, but the trend in 
p{T) remains similar to h/t = 0. Between h/t = 0.01 
and h/t = 0.02, however, there is a drastic change in 
p(T), and, as we will see later, in the magnetic state. 
The primary effect is a sharp increase in the T < T^ 
resistivity, with the T — > resistivity now being almost 
40% of the paramagnetic value. Even at this stage it is 
clear that p{T, h)/p{T, 0) can be very large as T — ^ 
and is -- 4 for T - T° and h/t = 0.02. Increasing the 
field even further leads to a reduction in p(T) over most 
of the temperature window since the field promotes a 
ferromagnetic state suppressing the AFM fluctuations. 

We can study the impact of the magnetic field also 
within the ZFC scheme. Fig. 2 shows the result of apply- 
ing a field after the cooling the system to four different 
temperatures, (i) T+ (shghtly above T°), (n) T" (slightly 
below T^), (in) Tf/2 and (iv) T°/A. 

For T > T^, the applied field mainly suppresses the 
AFM fluctuations, leading to a gradual fall in the resis- 
tivity. This is weak negative MR. Below T" (where there 







FIG. 3: Colour online: Evolution of spin correlations in the 
clean system in response to a magnetic field. The plot shows 
fi = So. Si in magnetic snapshots at different h, where So is a 
reference spin (bottom left corner). T — T^ /2 and the fields 
are, from left to right, h/t = 0, 0.008, 0.02, 0.05. 



is already noticeable AFM order) and at T"/2 the resis- 
tivity remains almost unchanged till some value hc{T) 
and then switches dramatically. There is a peak in the 
ratio p{T,h)/p(T,0) at h > 0.02 and then a gradual de- 
crease. At T'°/2 the ratio reaches a maximum ~ 12. 
At lower temperature, T°/4, the field appears to have a 
much weaker effect, mainly because the update mecha- 
nism that we adopt does not allow a cooperative switch- 
ing of the AFM state at low T till very large fields. Over- 
all, there is a window of T over which a moderate mag- 
netic field can lead to a several fold rise in resistivity. 

A first understanding of the rise in resistivity can be 
obtained from the magnetic snapshots of the system at 
T = T^/2 in Fig.3. We plot the correlation fi = Sq.S, 
in an equilibrium magnetic snapshot, where Sg is a ref- 
erence spin (bottom left corner) and Sj is the spin at 
site Hi. The left panel is at h/t = and shows a high 
level of A type correlation (stripe like pattern). This 
is a 'low dimensional' electron system since the electron 
propagation is along the one dimensional stripes in this 
2D system. The stripe pattern has a high degree of or- 
der so the scattering effects and resistivity are low. The 
second panel is at h/t — 0.008, just below field induced 
destruction of AFM order, and the pattern is virtually 
indistinguishable from that in the first panel. 

The third panel in Fig.3 is at h/t = 0.02 where the 
applied field has suppressed long range A type order. 
However, there are strong A type fluctuations that per- 
sist in the system and they lead to a pattern of short 
range ordered A type patches with competing orienta- 
tions, {Qai;Qa2} and {Qa3:Qa4}j in a spin polarised 
background. This patchwork leads to a high resistivity, 
higher than that in leftmost panel, since the ferromag- 
netic paths are fragmented by intervening A type regions, 
while the A type regions are poorly conducting due to 
their opposite handedness. In the last panel the field, 
h/t = 0.05, is large enough so that even the AFM fiuc- 
tuations are wiped out and the spin background is a 2D 
ferromagnet with extremely short range inhomogeneties. 
The resistance here is significantly below the peak value. 

Let us summarise the physical picture that emerges in 
the non disordered system before analysing the effect of 
antisite disorder. The ingredients of the large MR are the 



following: (i) An AFM metallic phase, without too much 
quenched disorder so that the resistivity in the magnet- 
ically ordered state is small, (ii) Field induced suppres- 
sion of the AFM order at h = hc{T), say, replacing the 
ordered state with AFM correlated spins in a FM back- 
ground, leading to a high resistivity state. In contrast 
to the standard negative MR scenario, where an applied 
field pushes the system from a spin disordered state to 
a spin ordered state, here the applied field pushes the 
ordered (AFM) state towards spin disorder. While our 
concrete results are in the case of a 2D 'double perovskite' 
model and a stripe-like ground state, the principle above 
is far more general and should apply to other non ferro- 
magnetic ordered states in two or three dimensions and 
to microscopic models that are very different from the 
double perovskites. We will discuss this issue at the end. 

Defects are inevitable in any system and in particular 
one expects antisite disorder in the DP's. The concen- 
tration of such defects may actually increase on electron 
doping a material like SFMO, due to the valence change, 
and we need to check if the large MR is wiped out by weak 
disorder. The presence of ASD affects the zero field mag- 




FIG. 4: Colour online: Field response in the presence of an- 
tisite disorder. The temperature is T = T'°/2, where T^ is 
the Tc at h/t = in the non disordered system. The struc- 
tural order parameters are indicated in panel (c) , S = 1 — 2a:: 
where x is the fraction of mislocated sites. The legend is same 
for all the panels. The results are obtained via ZFC. Panel 

(a) Field dependence of resistivity, normalised to p at h/t — 0. 

(b) Magnetic structure factor at the major AFM peak Qai. 
The value is same at Qa2 also, (c) Growth in the FM struc- 
ture factor with h. (d) Growth in the complementary AFM 
peak Qa3, the result is same for Qa4. The primary observa- 
tion here is the continuation of the clean limit results, with 
suppressed magnitude, down to 5* = 0.50. 




FIG. 5: Colour online: Field response in the antisite disor- 
dered systems at T = 7° /2, for the S values in Fig. 4. The 
left panels indicate the structural domains in the antisite dis- 
ordered systems. The middle column shows the spin correla- 
tions at h/t = 0, note that rows 1-3 show significant A type 
order, while the pattern in the 4th row has AFM domains of 
both orientations and so suppressed long range order. The 
right column shows the spin correlations at h/t — 0.02. In 
rows 1-3 the AFM pattern gets fragmented and FM regions 
show up. In row 4 the finite field pattern is not significantly 
different from the h/t = case. Overall, the field enhance- 
ment of 'spin disorder' is large in the first three cases but 
modest at strong ASD. 



netic state itself, as we have discussed elsewhere [18] , and 
the field response has to be understood with reference to 
this h/t = state. Let us first define the augmented 
model to describe the presence of ASD. 

We generalise the clean model by including a config- 
uration variable ry^, with rji = 1 for B sites and rji = 
for B' sites. In the structurally ordered DP the rji al- 
ternate along each axis. We consider progressively 'dis- 
ordered' configurations, generated through an annealing 
process [2U] . These mimic the structural domain pattern 
observed [211 122] in the real DP. For any specified {rj} 
background the electronic-magnetic model has the form: 

H == Hioc{l]} + Hkin{r/} + Hmagiv} 

Hiociv} = es Ei VifLfio- + <^B' J2i (1 - Vi) ml^m.^, the 
hopping term connects NN sites, irrespective of whether 
they are B or B': Hktn{v} = -tiY.{r,j)^ViVjfLfj'y 



'^Viilj){fia''^j<^ + h.c). The magnetic interactions include 
the Hund's coupling on B sites, and AFM superexchange 
Jaf — 0-08t between two NN B sites: Hmagiv} — 

JJ2iVi^i-fL^ai3fii3 + JAFj2{i.j) ViVj^f^j- For simplic- 
ity we set the NN hopping amplitudes ti = t2 = is = t- 

We use the simplest indicator 5 = 1 — 2a; (not to be 
confused with the magnitude of the core spin), to char- 
acterise these configurations, where x is the fraction of 
B (or B') atoms that are on the wrong sublattice. Fig. 4 
shows the resistivity ratio p{h)/p{0) at T = T"/2, in 
panel (a), and the field dependence of structure factor 
peaks in panels (b)-(d). Fig. 5 first column shows the 
structural motifs on which the magnetism is studied. 

Down to S" = 0.50 the ratio p{h)/p{0) has a behaviour 
similar to the clean case. Fig. 2, but the peak ratio re- 
duces to ^ 3 for S — 0.50. There is a corresponding 
field induced suppression in the principal AFM peak Qai 
and an enhancement of the FM peak Qp. The comple- 
mentary AF peak Qas slowly increases with h, has a 
maximum around h/t = 0.02 (where the disconnected 
AFM domains exist) and falls at large h as the system 
becomes FM overall. The trend that we had observed in 
the clean limit is seen to survive to significant disorder. 
At S" = 0.08, where the B-B' order is virtually destroyed, 
the h/t = state. Fig. 5 last row, has no long range AFM 
order. It is already a high resistivity state and an applied 
field actually leads to weak negative MR. 

Let us place our results in the general context of AFM 
metals, (i) Theory: we are aware of one earlier effort 
[23] in calculating the MR of AFM metals (and semi- 
conductors), assuming electrons weakly coupled to an 
independently ordering local moment system. Indeed, 
the authors suggested that AFM semiconductors could 
show positive MR. Our framework focuses on field in- 
duced suppression of long range AFM order, rather than 
perturbative modification, and the positive MR shows up 
even in a 'high density' electron system. The electron- 
spin coupling is also (very) large, J/t 3> 1, and cannot 
be handled within Born scattering. 

(ii) Experiments: the intense activity on oxides has led 
to discovery of a few AFM metals, e.g, in the manganite 
Lao.46Sro.54Mn03 [5], in the ruthenate Ca3Ru207 [6l|7], 
and CaCrOa [8]. Of these Ca3Ru207 indeed shows large 
increase in resistivity with field induced growth of FM 
order [7]. For the manganite and CaCrOs (where the 
resistivity is too large) we are not aware of MR results 
across the field driven transition. Some heavy fermion 
metals too have AFM character but the field driven tran- 
sition affects the local moment itself, making the present 
scenario inapplicable. 

(iii) Generalisation: the information about the spin 
configurations at any (T, h) is encoded in the structure 
factor, -D(q). Let us put down a form for D{q) and 
suggest how it affects the resistivity. To simplify notation 



we will assume that the AFM peak is at one wavevector 
Q, while the FM peak is at {0,0}. The AFM phase 
has an order parameter itiaf, say, while, for h > hc{T), 
there is induced FM order of magnitude mp . Assuming 
that the dominant fluctuations in the relevant part of the 
(T, h) phase diagram are at q ^ Q, we can write: D{q) ^ 
m\p5{q_-q,) + A/{l + {q_-q,f(^) when /i< K{T), and 
£)(q) ~ ml6{q) + A' / {I + (q - Q)^^^) ^Yi^^ ^ > ^^i^j^y 
mAFi'mF and ^ depend on {h, T), A depends on rriAp and 
^ and A' on mp and ^. The amplitudes A and A' vanish 
as the corresponding m tend to saturate (since the spins 
get perfectly ordered). The delta functions D(q) dictate 
the bandstructure while electron scattering is controlled 
by the Lorentzian part. 

Consider three cases (a) T = T+ , h = 0, {h) T = 
T-, h = 0, and (c) T = T- , h> hc{T). In (a) there is 
no order, so A is large, and p ~ pa, say. For (b) even if 
^ were the same as in (a), the presence of a large order 
parameter would suppress A and hence the scattering. 
We call this p = p^ « Pa (assuming there is indeed 
a large order parameter and no significant background 
resistivity due to impurities) . If we apply a field such that 
niAF — > but mp is still small, then the structure factor 
crudely mimics the paramagnetic case, and we should 
have pc ~ Pa- If all this is true, then just beyond field 
suppression of AFM order (and for T just below T^) we 
should get Pel Pb ^ 1- Broadly, if the appearance of 
AFM order with reducing T leads to a sharp drop in p 
then the field induced resistivity ratio can be large. This 
is independent of dimensionality and microscopic detail. 
A caution: as T — >■ 0, the applied field would drive a 
first order transition from a large m,AF state to one with 
large mp and weak scattering. The scenario above will 
not work, as our Fig. 2 illustrated. 

Conclusion: We have studied the magnetoresistance 
in an antifcrromagnetic metal motivated by the predic- 
tion of such a phase in the double perovskites. Beyond 
the modest field needed for suppression of long range an- 
tifcrromagnetic order, the system shows almost tenfold 
increase in resistivity near Tc- The effect originates from 
strong antifcrromagnetic fluctuations in the field induced 



ferromagnetic background. The large positive magne- 
toresistance, though suppressed gradually, survives the 
presence of significant antisite disorder. The principle 
that we uncover behind this 'colossal positive magne- 
toresistance' should be applicable to other local moment 
based AFM metals as well. 
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